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Abstract 

A generalized algebra of noncommutative coordinates and momenta embracing non-Abelian 
gauge fields is proposed. Through a two-dimensional realization of this algebra for a gauge field in- 
cluding electromagnetic vector potential and two spin-orbit-like coupling terms, a Dirac-like Hamil- 
tonian in noncommutative coordinates is introduced. We established the corresponding energy 
spectrum and from that we derived the relation between the energy level quantum number and 
the magnetic field at the maxima of Shubnikov-de Haas oscillations. By tuning the noncommuta- 
tivity parameter 9 in terms of the values of magnetic field at the maxima of Shubnikov-dc Haas 
oscillations, we accomplished the experimentally observed Landau plot of the peaks for graphene. 
Accepting that the experimentally observed behavior is due to the confinement of carriers, we con- 
clude that our method of introducing noncommutative coordinates provides another formulation 
of the confined massless Dirac fermions in graphene. 
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1 Introduction 

The recent experimental observations of the anomalous quantum Hall effect in monocrystalline graphite 
films of one atomic layer thickness [H [2] revealed the fact that in this material, called graphene, elec- 
trons behave as effectively massless relativistic particles. Theoretically, this unexpected quantization 
of Hall conductivity can be explained in terms of the massless Dirac-like theory O H]. On the other 
hand magnetic oscillations of electrical conductivity known as Shubnikov-de Haas (SdH) oscillations [5] 
were measured in patterned epitaxial graphene 0. It was shown that its transport properties result 
from carrier confinement and coherence. Moreover, in [6], to explain the observed behavior of the 
maxima of SdH oscillations, an analytic expression for the energy levels which takes into account the 
confinement of the charge carriers due to the micrometer-scale of the sample, has been proposed. In 
fact, it is in accord with the theoretical study of confining massless Dirac fermions by introducing a 
coordinate dependent mass term [7]. 

Noncommutativity of coordinates naturally appear in Landau problem as well as in string theories 
where noncommutativity is proportional to the magnetic field in the former and to the background 
fields in the latter. However, one can also introduce noncommutativity of coordinates by the con- 
stant noncommutativity parameter 9 as an intrinsic property of the space. Then it is legitimate to 
consider noncommutative version of any physical system employing appropriate formulation. When 
noncommutativity of space is imposed without attributing a definite physical meaning to the noncom- 
mutativity parameter 6 one should provide its interpretation. Obviously, noncommutativity can be 
taken as a fundamental property of space by claiming that the noncommutativity parameter 9 is very 
small and it is responsible for the errors in the measurements of the related physical quantities. In 
this case each physical system will provide a different limit for the possible value of the noncommuta- 
tivity parameter 9 (see e.g. |8] and references therein). However, there exists another interpretation: 
Noncommutative formulation of a dynamical system can be considered as a tool to link some diverse 
manifestations of a physical phenomena, e.g. obtaining the fractional Hall effect from the Hall effect 
in noncommutative coordinates by an appropriate choice of the noncommutativity parameter 9, as 
was reported in [9]. Because of retaining the terms up to a definite order in 9 one should show that 
in the dynamical problem considered the fixed value of 9 is in accord with this approximation. This 
interpretation of noncommutating coordinates as a linkage between different phenomena may give 
same clues of finding an easier method of formulating interacting systems from noninteracting theo- 
ries [To]. Hence, we would like to study massless Dirac theory in two-dimensional noncommutative 
space to perceive whether a similar interpretation of noncommutativity is possible which can yield a 
better understanding of some peculiar properties of graphene. However, introducing noncommutative 
coordinates into spin-dependent Hamiltonian systems defined by constant matrices is not well estab- 
lished. The usual method of introducing noncommutativity is to replace ordinary products with star 
products which is equivalent to the shift 



where are the quantum phase space variables. Obviously, this method, which does not take 

into account spin degrees of freedom, is not suitable to deal with matrix valued, constant observables. 
Utilizing the semiclassical techniques developed in a method of introducing noncommutativity 
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of space appropriate to deal with non-Abelian spin matrices was presented in |12[ [5], However, this 
semiclassical treatment is not amenable to introduce noncommutativity into Dirac equation of fermions 
interacting with external fields. We propose a solution for this issue. First quantum commutation re- 
lations between phase space variables in noncommutative space suitable to consider non-Abelian fields 
are presented. We derived them from the semiclassical brackets of classical phase space coordinates 
proposed in [12]. Then, realizations of this algebra may be employed to introduce noncommutative 
coordinates in Hamiltonian systems. It is applicable to systems where the interaction terms are co- 
ordinate independent but non-Abelian because of being matrices. Obviously, this formalism can be 
applied to Hamiltonian systems whose interaction terms are not matrices, so that it establishes an 
alternative to the custom star product approach which is equivalent to the shift ([T|). We then focus 
on the confinement problem of massless Dirac particles in graphene and propose that introducing 
noncommutative coordinates within our approach may be used to deal with one of its basic features: 
the Landau plot of the maxima of SdH oscillations. To obtain the 6 deformation we give a realization 
of the generalized canonical commutation relations in the presence of spin-orbit-like couplings and a 
transverse magnetic field. We obtain the spectrum of the proposed Hamiltonian and show that it can 
be used to formulate the SdH effect in graphene. By fixing the noncommutativity parameter 9 we 
actually established a good agreement with the experimental observations which are known to result 
from the confinement of massless Dirac fermions. In the last section we discuss the obtained results 
and their future applications. 



2 Generalized algebra 

We would like to present a formulation of quantum mechanics in noncommutative coordinates acquired 
by quantizing the semiclassical approach of [11]. Hence, let us briefly recall this formulation. When one 
considers a dynamical system in the classical phase space (P/, Qi); I = 1, ■ ■ ■ , M one can introduce 
the star product 



exp 



ih 



(2) 



2 yOQ^ dPi dPidQ^ 

to achieve quantization of the system within the Weyl-Wigner-Groenewold-Moyal approach [13]. Cor- 
responding to the quantum commutators one considers the Moyal bracket of the observables f{P, Q) 
and g{P-, Q) which are some functions as 

[/(P, Q),g{P, Q)], = /(P, Q) * g(P, Q) - g{P, Q) * f{P, Q). 

Poisson brackets follow in the classical limit, 

[/e-.W.ff^-Wl. = UiP^QUiPm ^ - §;^. (3) 

For the matrix observables Mki{P, Q) and Nki{P, Q) the Moyal bracket can be defined as 

([M(P, Q), iV(P, Q)] J,, = Mkm{P. Q) * N„a{P, Q) - Nkrn{P, Q) * MrruiP, Q)- (4) 

In contrary to the observables which are functions, matrix observables in general are not Abelian, 
so that the classical limit ^ of their Moyal bracket (HD will yield, in general, a term behaving as 
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which would be singular. However, if the observables are related to spin they may possess H 
dependence, then the singularity which we mention is superfluous. To take into account this fact we 
define the "semiclassical" bracket 

^[M,N] + ^{M{P,Q),N{P,Q)} - J 



{M{P, Q), NiP, Q)}c = -^[M, N] + -{M{P, Q),N{P, Q)} - -{N{P, Q),M{P, Q)}, (5) 



where one should retain the terms up to h. The first term on the right hand side is the ordinary 
commutator of the matrices and the last two terms are Poisson brackets. The bracket ([5]) does not 
satisfy Jacobi identities. This is due to the fact that in its definition one keeps the first two terms of 
the fully-fledged Moyal bracket (jl]). However, in our approach the semiclassical limit is taken after 
multiplying the observables by the star product ([2]). Hence, we should consider the semiclassical limit 
of the Jacobi identity which is the first two terms of the Moyal bracket relation 

-■^({M,{iV,L},}, + {iV,{L,A/},}. + {L,{i\/,iV}J,) = 0_i + Oo{h'') + Oi{h) + ■ ■ ■ . 

In fact, one can show that the semiclassical limit of the Jacobi identity is satisfied 

O_i(i) + Oo(;i°) = -■^[M,[iV,L]] + [M,{iV,L}]-[M,{L,iV}] + {M,[iV,L]} 

-{[iV, L], M} + (cyclic permutations of M, iV, L) = 0. 

Nevertheless, once we perform quantization and deal with quantum operators choosing a realization 
of quantum phase space variables, we should impose that they satisfy Jacobi identities. 

The first order matrix Lagrangian adequate to formulate spin dynamics in noncommutative coor- 
dinates is 



2 - ■ - . HQ{r,p) (6) 

where a, /3 = 1, • • • , d. We would like to emphasize that Aa is in general matrix valued, p denotes the 
related coupling constant and I is the unit matrix. The constant, antisymmetric non-commutativity 
parameter Oap appears divided by h to set its dimension at (length)^. The definition of canonical 
momenta 

Vr" = — , TT" = — 



yields the dynamical constraints 



V''" = (7) 



= U- + -r- I+^p^. (8) 



By setting P/ = (vr", vr^ ), Qi = {r°',p°^) in ([5]) one can show that the constraints (I7| and ([8]) obey the 
semiclassical brackets 

{V'a,^/3}C = pFap, 
r / 2 ;2i 



h 



ah 
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5ai3 is the Kronecker delta and Fq,^ is the field strength, 

dAfs dAa, ip 

where the last term is the ordinary matrix commutator. Thus, we may classify ip^\ z = 1,2, as second 
class constraints and the matrix whose elements are 

Q' = {C,V'|'}c, (10) 

possesses the inverse C~^: 

C^J^'C-Jf = 5i5l,. (11) 
The inverse matrix elements can be employed to define the "semiclassical Dirac bracket" as 

{M, N}cD = {M, N]c - {M, r]c C-}, {V'"', iV}c, (12) 

so that the constraints ([7]) and dH) effectively vanish. The basic classical relations between the phase 
space variables following from ([6]) can be established, at the first order in Q and keeping at most the 
second order terms in p, as 

K,r/^}CD = (13) 

{p",/}cD = pF'^f'-^iFerr^, (14) 

K,/}cz) = 6"^-^{eFr^ (15) 

where (flF)"^ = O'^^F^, {eF9)°'^ = 9°'^F^e^. We omitted the identity matrix I on the left hand sides. 
Indeed, in the sequel we will not write I explicitly. 

The brackets (|13p -(|15 p differ from the Poisson brackets up to commutators of matrices, so that 
for observables which are not matrices they reduce to the ordinary Dirac brackets. Therefore, we can 
extend the canonical quantization rules to embrace the matrix observables by substituting the basic 
brackets with the quantum commutators as { , }cd ~^ ^[ > ]g- To distinguish the matrix commutators 
and quantum commutation relations we denoted the latter as [ , ]q. This yields the generalized algebra 

K,f% = i^^ (16) 
= ihpF'"^ -ip'^{F9F)'^^, (17) 

[f",/]q = m'^^ -ip{eFY^, (18) 

[p°,r/3]g = -ins'^P +ip^FeY^. (19) 

Note that, on the right hand side we keep the first order 6 contributions, so that everything can 
only depend on x^, defined as r^lg^o = a^a- For Abelian gauge fields this type of algebra has already 
been considered in [l3| and a similar one in noncommutative space for an electromagnetic field was 
discussed in [15] (see also [16] and references therein). 

One can employ realizations of the algebra (|16p-()19p to introduce noncommutative coordinates. 
To illustrate it, let us deal with the commutative case ^ = and let the gauge field be not a matrix 
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but the 2-dimensional electromagnetic one a, = {—Br2/2, Bri/2), which leads to a constant magnetic 
field transverse to the (ri, r2)-plane B. For these choices the algebra becomes 



[^i,?^i]g = 0, \pi,pj\q = iehBeij, [fi,pj]g = ih5ij. (20) 
A realization of the algebra ()20p is 

d 

Pi = + ri = ri. (21) 

Through the substitution of classical momenta with the realization (j21|) . in the free Hamiltonian 
Hq = p^/2m, the minimal coupling to the gauge field in quantum mechanics can be achieved as 

1 / 3 \ ^ 
Hi^t = Ho{p,q) = — (-ih—+ea] . (22) 
2m \ or J 

We will extend this point of view to define quantum mechanics in noncommuting coordinates. 

Although we will employ another realization to propose a Hamiltonian adequate to describe 
graphene on the noncommutative plane, let us present a realization of (jl6P "(|19 p . In terms of the 
covariant derivative 

d 

Da = -ih- pAa = -ihVa - P^a, (23) 

OXa 

we can realize the algebra (fT6l) - (fT9]l by setting 

Pa = Da-^Fa^d[S^D^, (24) 

r„ = Xa- -^OapDp, (25) 
as far as F^p are constant, commuting matrices which are equivalent to the conditions 

- ihVaFp^ - p[Aa,Fp^] = 0, [Fap, F^s] = 0. (26) 



These conditions are also necessary to show that the realization (j24p -" (j25p satisfies the Jacobi identities. 
We would like to emphasize that this realization is valid for either Abelian or non-Abelian A^- It can 
be employed to introduce the related dynamical system in noncommutative coordinates as 

Hie)^H^o)ir,p), (27) 

where if(o) ii"^ p) is the free Hamiltonian appropriate to the considered system. Indeed, this constitutes 
an alternative method to the star-product approach of introducing noncommutative coordinates in 
quantum systems. 

3 Dirac particles in noncommutative space 

In graphene, around each Dirac point, which is the point at the corners of Brillouin zone, the free 
Hamiltonian is written as the massless Dirac-like Hamiltonian [17\ [T8] 

Hg\p,q)=Vfp-a (28) 
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for low energies and long wavelengths. Here, p = {px,Py) is the two-dimensional momentum operator 
and a = (cr^,cry) where crx,y,z are the Pauli matrices acting on the states of two sublattices. Vf is the 
Fermi velocity playing the role of the speed of light in vacuum. 

We would like to deal with the dynamics of the massless Dirac particle on the noncommutative 
(x, y)-plane whose free Hamiltonian is (j28|) within the method presented in the previous section. For 
this purpose let the gauge field be 



1,2, 



(29) 



which is non-Abelian. The first term corresponds to the transversal, constant magnetic field B and 
the others are spin-orbit-like coupling terms. However one should keep in mind that for graphene a 
act on the states of sublattices, so that though k and I, respectively, like the coupling constants related 
to the Rashba and Dresslhauss spin-orbit interaction terms for electrons, their effect is to give rise to 
terms proportional to o"^ and unity in the non-deformed Hamiltonian. In fact, shifting momenta in 
(j28p with the gauge field (j29p yields the following Dirac-like Hamiltonian 



Hd = Vf:{p-A)-a 

= v^ai i-ihVi + ^eijXj ) - 2vf (ka;, + /). 



(30) 



We would like to get a ^-deformation of this Hamiltonian employing the procedure outlined in the 
previous section. Hence, we set p = 1 and by using the definition ([9]) we obtain the field strength 
corresponding to (f29]) as 

(31) 



The algebra (jl6p - (|18|) now becomes 



ih[eB + - k^)a, ) e^j + ( ie^B'^O + ^eBe{l'' - k')a. 



Ai 



1 + + (/2 _ k2y_^_^ 



h 



B 



20 



(32) 
(33) 

(34) 



B — We deal with small /, k, so that we neglect the terms at the order I'^k"^ for n + m > 4. 



where /| 

Obviously, ([29]) and (f3T|) do not satisfy the conditions (126]) . so that one cannot make use of the 
realization (|24p . (|25p . Nevertheless, we accomplish a realization of (|32 p - (|34p as follows: 



Pi 



1 + i + 



21 



B 



eB 

-ihVi H — —^ijXj — ikeijGj — la^ 



26 



+(/ -k )-^enmXn{-ihS/m){ikeijaj +lai), 



1 + 



21 



2 

B 



Xi 



2h"'' 

„2 



eB_ 



(35) 



(36) 



One can demonstrate that (j35]) and ([36]) satisfy the Jacobi identities at the first order in 9 and ignoring 
the terms at the order of l^-k^ for n + m > 4. 
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Through the procedure outhned with (j27|) ^-deformation of the Hamiltonian (j3U|) can be achieved 
by substituting momenta with the realization (j35|) in the Hamiltonian (|28p as 



Trie) ^ 
° 2 



p- a + [p- a 



(37) 



Indeed, plugging (f35]l into (f37|l yields 



Vf{l + 



2/2 



-ihVi + 



e5 



ejj2;j I cj,: - 2t!F 



2Z 



B 







+ 2{l - k )^€nmXni-i^^m) {kaz + l). 



(38) 



One can observe that for ^ = 0, it yields the Hamiltonian given in (|30p . In [19j noncommutative 
structure emerging in graphene was studied where the noncommutativity parameter considered is due 
to a lattice distortion term present in the Hamiltonian. In fact, it is similar to the mass term, so 
that it is related the our constant parameter k. However, it is like a coordinate dependent mass term 
because they also consider a nonconstant distortion. 

To establish energy eigenvalues it is convenient to write (|38p in terms of the complex variables 
z = X + iy, z = X — iy as 



(0) 

D 



5'^ 



h 



iK [2hVz + ^z) 



iK {-2hV^ + ^z) 



(39) 



where = —iheijXiV j = h{zV ^ — zVz) is the angular momentum operator. The involved constants 
are defined as 



g± = -2v^{l±k) 



211 ~ 



h±=AvFj^{l±k){l' -k'), K = ^;f(1 + _) 

B 



To derive the eigenvalues of (j39p algebraically, we introduce two pairs of annihilation and creation 
operators: 



a 
b 



eB ■ 



-||:(2W, + f.), at = ^(-2W,+ 



which are mutually commuting and satisfy the commutation relations 

[a,at] = [6,6t] = i. 

Hence, the Hamiltonian (j39p acquires the form 



h+{b^b-a^a) Ka) 

ka g_ - h_{b^- a^a) 



where K = 2vf heB ( 1 + 



2ii 



The eigenvalue equation for the two component spinor 



(9) ( V'l 



E 



1p2 
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leads to two coupled equations 

g+ - h+{b^b-a^a) - E 



9' 



h^{b'^b- al^a) - E 



V'2 



-Kaipi. 



(40) 
(41) 



After some calculation, one can show that the equation satisfied by the spinor component ipi takes 
the form 



E^ + AE {kI + ^/(/2 - k^){l + 26+6 + 2a^a)^ - K^a^a + iK^if - k^) 



^1 = 0. (42) 



To draw the energy eigenvalues, let us write the state corresponding to the spinor component tpi as 

1 



=1 n, m) 



:(6t)"^+"(at)" I 0); 



(43) 



y^nl{m + n)! 

n, m = 0, 1, 2 • • • , and by definition a|0 >= 6|0 >= 0. In the complex plane ()43p yields 

{z,z I n,m) = Nmnz'^Ll' (y) e-^^^" 

where are the Laguerre polynomials and Nmn are the normalization constants whose explicit forms 
are not needed. 

Obviously, (|43p satisfies the relations 

(6^6 — a^a) \ n, m) = m \ n, m), 
a^a \ n,m) = n \ n,m) , 

where m and n are the quantum numbers corresponding, respectively, to the angular momentum 
eigenvalues and the Landau levels. Now, (|42|) can be solved to deduce the energy spectrum as 



En,Uk,l,0,B) = ±2vf (1 + ^ 



21 



B 



21% 



n 



(44) 



Moreover, one can show that the corresponding spinor components are given by 
with the convention V'-i.m = 0. Here s' is a constant, which can be read from (|4ip . 



I n, m) 
s' I n — 1, m + 1) 



4 Shubnikov-de Haas effect 

The Shubnikov-de Haas (SdH) effect is a magnetotransport phenomena that occurs in materials in 
a strong magnetic field of about 1 Tesla and for low temperature about few kelvins [S]. It is an 
oscillatory dependence of the electrical resistivity of a metal or a semiconductor as a function of the 
applied constant magnetic field. More precisely, the SdH effect is produced by the oscillations of the 
density of states at the Fermi level. The mechanism can be understood for metals considering Landau 
levels |2Uj which are the energy levels of electrons in the presence of a magnetic field. If the electrons 
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fill the energy levels up to the level n + 1, the Fermi energy which is equal to the chemical potential 
at absolute zero, will lie in this level. As the magnetic field increases the degeneracy of the Landau 
level increases. Thus the electrons move to the level n, depopulating the level n + 1, so that the 
Fermi energy is decreased. Now, increasing magnetic field leads to less populated Fermi level until 
all electrons migrate to the lower energy level. Hence, the conductance or the resistivity will oscillate 
as a function of the external magnetic field. The maxima of the SdH effect occur at the magnetic 
fields B]\f which can be calculated by equating the energy level corresponding to the index with 
the chemical potential fi (Fermi energy). Hence, the relation between N and -Bat predicted by our 
approach is established as 



N ' 



2ehBN 



4 + 4ii/ + 4(/2-A;2) 



-/(A:^-/2)(2?n + l) + -^+/ 



(45) 



by solving the equation Ei\f^rn{k,l,0, B^) = fi obtained from ()44p . For convenience, we rescaled the 
noncommutativity parameter as 

e = -'-^e. 

To analyze the SdH effect in graphene within our formulation we shall choose the involved param- 
eters adequately. To start with, we require that the spin-orbit-like coupling constants obey 

l = -J^ + k. 

With this choice ()45p is simplified and takes the form 

e /B{m,k) 



where we defined 

Since the noncommutativity parameter is a free parameter, it can be fixed in diverse fashions. 
However, one should keep in mind that its value should be consistent with the approximation of 
retaining the terms up to the first order in 6. In particular, for the limiting values of Sat, we propose 
to choose 9 as 

1 iBn, BN<^B{m,k)/k J 



where 7, /3 are two constants and we assume that k ^ 0. We can analyze (|46p separately for each case 
given in (j47p . For i?jv > B{m, k)/k we deduce the behavior 

N> = (48) 
n"' Bj^ 

by neglecting a term behaving as 1/Bj^. Thus, for large Bpf, N changes linearly with respect to l/Bj^. 
However, in the second case, B^ <C B{m, k)/k, N leads to the constant value 

^B(m,k) , , 

iV< = ^ (49) 
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Let us link these considerations to the experimental observations of [B] . They obtained the limiting 
values 

= 1 25, « 2.5 T I 

where the constant is given by 



This fixes the ratio 



— w 34 X lO^^'^ kg.m/s. 



Now, we would like to determine the value of the noncommutativity parameter 6 comparing (j50p 
with (j48|) and ()49p for m = 0. The other values of m can be treated similarly. First of all observe that 
we may impose 

B(0,k) , , 

\' ^ = 2.5 T. 51) 

k 

To simplify let k = {^/2vf)6, so that (I5ip yields the equation 

25^ - ( 3 + ^ V + 1 = 

-DO/ 



whose solutions are 
Hence, we may set 
which implies to choose 



0.77 ±0.3. 



= 1 X 10"^^ kg.m/s 



/3 f« 4 X 10"^^JmsT, 7 « 1 x 10"^^JmsT"\ 

It worths to observe that the magnitude of the noncommutativity parameter for the limiting cases 
(07]) reads 

\e{B)\ = [ ^iv' X 10^'' i?^>2.5T 1 
\ X 10^16 j^2^ S7v<2.5Tj' 

Therefore, there is no conflict with keeping the terms up to the first order in 9. Until now we dealt 
with the values of 6 for the limiting values of the magnetic field B^. However, we can also choose it 
appropriately for all values of Bjy. Inserting the choice (|5ip into (|46p yields 

0k /2.5 T \ 
N=To^ + l]- (52) 



n? \ Bn 

To write the full expression for 6, let us introduce the Heaviside step function 



H{x) 



0, 

1/2, 
1, 
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Figure 1: Landau plot of the maxima of SdH oscillations. 



which can be given analytically as |21] 

H{x) 



lim 



2 vr 

We choose the noncommutativity parameter to be 



11 .X 
- + - tan ^ - 



1 + 2.5B 



N 



1 1 -1 
— I tan 

2 vr 



0.4 - B 



N 



0.01 



+ 



53 



0.9 + S 



N 



1 1 

— \ — tan 

2 vr 



0.83 - B 



-1 



N 



0.01 



+ 24.8 



1 1 

— I tan 

2 vr 



5 



-1 
TV 



1 1 

- H — tan 

2 vr 



0.83 



5 



AT 



0.4 



0.01 



0.01 



(53) 



which produces the limiting values correctly and with this choice (j52p yields Figure 1 (using MATH- 
EMATICA). Indeed, we have chosen ([53]) appropriately so that the Landau plot of the peaks given in 
Figure 1 matches well with the experimental one obtained in [B]. Moreover, one can check that the 
order of magnitude of the noncommutativity parameter is ~ 10~^^ m^, so that it is in accord with 
the approximation of ignoring the second order terms in 6. 



5 Discussions 



The results which we obtained are twofold: 

• A new method of introducing noncommutative coordinates into quantum mechanics is estab- 
lished. 

• An analytic method of obtaining confinement of massless Dirac particles in graphene is proposed. 

We introduced a generalized algebra of quantum phase space operators in noncommuting space on 
general grounds with momenta involving non-Abelian gauge fields. This constitutes an alternative to 
the custom method of introducing noncommuting coordinates by star products. It may lead to some 
new features of quantum mechanics in noncommutative coordinates. Moreover, it should be possible 
to extend it to field theory formulations. These are currently under inspection. 

We considered a two-dimensional space by a particular choice of gauge fields. A realization of 
the associated algebra is presented and employed to obtain a massless Dirac-like Hamiltonian on the 
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noncommutative plane. Its energy eigenvalues are established. Through an appropriate choice of the 
noncommutativity parameter 9 we showed that this energy spectrum is adequate to accomplish the 
experimentally observed behavior of the SdH oscillations in graphene, which are known to result due 
to the confinement of its charge carriers which are massless Dirac particles. Obviously, our main 
objective is to employ this noncommutative theory to understand those features of graphene which 
are not well understood within other formalisms. This work should be considered as the first step in 
this direction. We obtained a satisfactory noncommutative version of Dirac-like theory of graphene 
which led to some predictions. One of the next steps would be to obtain a field theory in terms of the 
Hamiltonian (I38p , which can be used to introduce other interactions like the spin of electron as in [22] 
into the noncommutative theory. 
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